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ABSTRACT 

An  important  unsolved  problem  in  nonlinear  control  is  that  of  stabilization  with  a  prescribed 
region  of  stability.  In  this  paper,  sufficient  conditions  are  obtained  for  the  existence  of  a  linear  feedback 
stabilizing  an  equilibrium  point  of  a  given  nonlinear  system  with  the  resulting  region  of  asymptotic 
stability  (RAS)  containing  a  ball  of  given  radius.  Conditions  for  global  stabilization  are  also  given. 
Feedback  stabilization  is  achieved  while  satisfying  a  certain  robustness  property.  The  technique  is 
applied  to  planar  systems,  resulting  in  a  complete  design  methodology  for  this  case.  Examples  and 
simulations  illustrating  the  method  are  presented. 

1.  INTRODUCTION 

Given  a  nonlinear  control  system  x  —  f(x,u )  where  /  :  IRn  x  IRm  — >  IRn  is  a  smooth  mapping  satisfying 
/(0,0)  =  0,  it  is  well  known  that  if  the  associated  linearized  system  is  stabilizable,  then  there  exists  a 
linear  feedback  control  which  also  stabilizes  the  (null  solution  of)  the  original  nonlinear  system.  The 
null  solution  of  the  closed-loop  nonlinear  system  of  course  possesses  a  nonempty  region  of  asymptotic 
stability  (RAS).  The  size  of  the  RAS  is  usually  not  stated  as  an  explicit  control  objective.  The  reason 
for  this  is  the  lack  of  systematic  analytical  tools  for  the  synthesis  of  feedback  control  laws  achieving 
specifications  on  the  RAS.  The  importance  of  obtaining  such  tools  is  clear,  and  has  been  emphasized 
in  [1]. 

The  traditional  approach  based  solely  on  linearization  at  an  operating  point  is  often  considered 
unreliable  from  a  stability  point  of  view.  An  alternative  design  method  consists  in  repeated  testing  of 
the  performance  of  the  closed-loop  system  for  each  of  a  set  of  possible  stabilizing  control  laws.  Since 
approximation  of  the  obtained  RAS  is  often  very  difficult,  each  of  these  tests  typically  involves  many 
simulations  of  the  closed-loop  dynamics,  and  the  method  is  hence  very  costly  [2]. 

In  this  paper,  sufficient  conditions  are  obtained  for  the  existence  of  a  linear  feedback  stabilizing 
an  equilibrium  point  of  a  given  nonlinear  system  with  the  resulting  region  of  asymptotic  stability 
(RAS)  containing  a  ball  of  given  radius.  Conditions  for  global  stabilization  are  also  given.  Feedback 
stabilization  is  achieved  while  satisfying  a  certain  robustness  property.  Synthesis  of  the  desired  feedback 
control  laws  rests  on  the  solution  of  certain  nonstandard  questions  in  linear  systems.  These  questions 
are  addressed  successfully  for  the  case  of  planar  systems,  for  which  a  complete  design  methodology  is 
achieved.  Examples  and  simulations  illustrating  the  method  are  presented. 


|  On  leave  from  Faculte  des  Sciences  de  Marrakech,  Universite  Cadi  Ayyad,  Marrakech,  Morocco. 


2.  SUFFICIENT  CONDITIONS  FOR  STABILIZATION 
WITH  PRESCRIBED  RAS 


Before  presenting  the  main  results  of  this  section,  we  establish  notation.  With  S  a  subset  of  (V,  ^Re(S) 
denotes  the  set  {92e(s)  :  s  £  S}.  For  a  real  matrix  M ,  [M],  and  [A/],,  denote  its  symmetric  and 
skew-symmetric  parts,  respectively: 

Ms  :=  1(M  +  Mt), 

Ms, 


The  spectrum  of  M  is  denoted  by  <r(M).  For  x  a  vector  in  7/2",  |x|  denotes  its  Euclidean  norm. 

0 

Denote  by  B(R)  the  open  ball  in  IRn  of  radius  R  centered  at  the  origin.  By  7/2_  and  (F_,  we  intend 
(—oo,0)  and  the  open  left-half  of  the  complex  plane,  respectively.  For  S  a  given  set,  Sn  denotes  the 
Cartesian  product  S  x  S  x  •  •  •  x  5  (n  times). 

Two  definitions  relating  to  stabilizability  of  linear  systems  are  now  introduced.  Let  A  and  B  be 
real  matrices  of  dimensions  n  x  n  and  n  x  m,  respectively. 

Definition  1:  Say  that  the  pair  (A,B)  is  symmetrically  stabilizable  1  if  3 K  €  IRmxn  such  that 
<t([A  +  BA],)  C  IRn_  •  For  A  a  nonempty  subset  of  2/2" ,  the  pair  (A,B)  is  said  to  be  symmetrically 
stabilizable  within  A  if  VA  €  A,  3 K  €  IRmxn  such  that  <r{[A  +  BI< ],)  =  A. 

Definition  2:  Sav  that  the  pair  (A,B)  is  normally  stabilizable  if  3 K  £  IRmxn  such  that  a(A+BK)  C 

On  On 

<F_  with  A  +  BK  a  normal  matrix.  Let  A  C  <F_  be  nonempty.  Say  that  (A, B)  is  normally  stabilizable 
within  A  if  VA  £  A,  3 K  £  IRmxn  such  that  <r(A  +  BK)  =  A  with  A  +  BK  normal. 


It  is  a  simple  exercise  to  show  that  if  ( A,B )  is  normally  stabilizable  within  A  then  it  is  symmet¬ 
rically  stabilizable  within  32e(A). 

We  consider  nonlinear  multi-input  control  systems  x  =  F(x)  +  Bu  where  F  :  IRn  — +  272"  is  analytic 
over  2/2"  and  satisfies  F(0)  =  0,  and  where  B  is  an  n  x  m  matrix.1 2  It  is  convenient  to  rewrite  the 
model  in  the  equivalent  form 

x  =  Ax  +  Bu  -f  h(x),  (1) 


where  A  := 
Let 


|A(q)  and  h(x)  represents  higher  order  terms. 


Ti  :=  {  h  :  IRn 


dh 


IRn ,  analytic  over  7/2"  :  /t(0)  =  0,  —■  (0) 

ox 


0} 


and  define  for  a  fixed  72  >  0  the  following  two  functions  on  H 
(i)  ||.||a  :  H  — ►  2/2+  such  that 


I  Mr  = 


sup 
*eB(R) 
Tjt  0 


IM*)I 

1*1 


1  Recall  that  Hurwitz  stability  of  [AT],  implies  that  of  M. 

2  Notice  that  this  model  is  not  restrictive  since  a  more  general  model  i  =  f(z,  v)  may  always  be  put 
in  the  form  above  by  letting  *  be  the  augmented  state  ( z,v)T  and  taking  u  =  ii. 


(ii)  ||.||  :  7i  —*  IR+  U  {00}  such  that 


*11  = 


sup 

T,alir 

x^O 


1M*)1 

1*1 


Notice  that  lini|x|_0  =  0  is  well  defined  and  that  we  allow  the  possibility  ||/i||  =  00  with  the 

understanding  that  inequalities  such  as  ||/ii  +  /isli  <  ||hi||  +  ||/i2||  are  to  be  interpreted  in  the  obvious 
way  when  ||/ii||  or  ||/i2||  is  infinite.  Then  (H,+, .)  with  either  ||.||fl  or  ||.||  is  a  normed  vector  space.  This 
can  be  seen  easily,  noting  that  analyticity  of  the  elements  of  %  implies  that  a  function  which  vanishes 
within  B(R )  must  also  vanish  everywhere.  We  naturally  define  two  balls  in  R: 


Br{p)  .=  {hen  :  ||h||/?  <  p }, 
B(p):={heH:  |H|  <  p}. 


Theorem  1:  Fix  R  >  0  and  let  ( A,B )  be  symmetrically  stabilizable  within  A  C  (—00,  —  ||/i||.r  )” 
(resp.  (—00,  —  ||fi||  )n).  Then  the  nonlinear  control  system  (1)  is  stabilizable  within  B(R )  (resp.  globally 
stabilizable)  using  linear  state  feedback. 

The  proof  of  this  theorem  relies  on  Proposition  1,  given  next.  Let 


x  =  F(x),  F(  0)  =  0  (2) 

where  F  is  analytic  over  IRn .  Let  the  null  solution  of  (2)  be  asymptotically  stable  in  the  sense  of 
Lyapunov.  Denote  the  associated  RAS  by  Vx .  Consider  a  change  of  coordinates  2  =  QTx  where  Q  is 
an  orthogonal  matrix  (i.e.,  Q  satisfies  QTQ  =  QQT  —  I).  In  the  new  coordinates, 

z  =  F(z)  (3) 

where 

F(z)  :=  QtF(Qz). 

Clearly,  the  origin  is  also  asymptotically  stable  for  Eq.  (3).  The  sets  Vx  and  Vz  are  in  general 
different.  However,  we  can  use  the  fact  that  orthogonal  transformations  preserve  norms  and  angles  to 
obtain  the  following  proposition. 

Proposition  1:  The  largest  Euclidean  balls  in  Vx  and  Vz  are  identical. 

Consequently,  for  each  R  >  0 


B(R)  CVX  <=>  B{R)CVZ. 


Proof  of  Theorem  1 :  Let  (Ai,  A2, . . . ,  A„)  £  A.  Since  (A,B)  is  symmetrically  stabilizable  within 
A,  there  is  a  feedback  gain  matrix  I\  such  that 

a([A  +  RA]3)  =  {A1,A2,...,An}. 

By  setting  u(x)  =  Kx  in  (1),  we  obtain  the  closed  loop  system 


x  =  (A  +  BK)x  +  h(x) 

=  [A+  BK}„x  +  [A  +  BK],x  +  h(x). 


(4) 


Since  [A  +  BK]S  is  symmetric,  it  can  be  diagonalized  using  an  orthogonal  transformation.  Let  Q  be 
such  a  transformation  and  define  new  coordinates  z  =  QT x.  Then  z  satisfies 


z  =  Gz  +  Dz  +  h(z )  (5) 

where 

D  -  diag(Aj,  X2, . . A„), 

G=Qt[A  +  BK]„Q 

and 

h(z)  =  QTh(Qz). 

Now  consider  the  Lyapunov  function  candidate  V(z)  =  L  |z|2  and  note  that  B(R)  is  a  level  set  for  V. 
The  derivative  of  V  along  trajectories  of  (5)  is 

V(z)  =  zTGz  +  ztDz  +  zTh(z). 

The  first  term  vanishes  since  G  is  skew-symmetric.  Thus, 

»=i 

<  max  (A,-)  \z\2  +  ||A||h  |z|2 

l<t<n 

=  (  maxn(A0  +  ||A||a  ^  \z\2  (6) 

for  all  z  €  B(R).  Noting  that  A  C  (— oo,  —  |J/r||/i)" ,  we  have  that  V(z )  <  0  for  all  nonzero  z  E  B{R). 
A  standard  Lyapunov  stability  result  now  implies  B{R)  C  Vz .  In  view  of  Proposition  1,  an  analogous 
statement  also  holds  for  Eq.  (4).  This  proves  the  first  assertion  of  Theorem  1.  The  second  assertion 
similarly  follows  from  the  observation 

v(z)  <  (  1I^f<xn(A«')  +  llftll  )  M2 

<0 

for  all  z  E  IRn ,  z  /  0.  This  proves  global  asymptotic  stability. 

□ 

Robustness  of  the  stabilization  property  of  Theorem  1  with  respect  to  perturbations  in  the  nonlinear 
terms  h(x)  is  now  considered.  The  higher  order  terms  h(x)  do  not  affect  asymptotic  stability  of  the 
null  solution  of  a  hyperbolic  system  (linearization  with  no  imaginary  eigenvalues).  In  our  framework, 
we  note  that  for  u(x)  =  I\ x  a  linearly  stabilizing  feedback,  the  null  solution  of  x  =  (A  +  BK)x  +  h(x)  is 
asymptotically  stable  for  all  h  EW.  However,  the  RAS  does  indeed  depend  on  variations  in  h.  The  next 
result  states  that  the  linear  feedback  u(x)  =  Kx  in  Theorem  1  is  robust  to  variations  in  h.  Specifically, 
the  assertion  is  that  B(R)  is  guaranteed  to  be  within  the  RAS  for  each  member  of  a  family  of  systems 
each  of  whose  linear  parts  is  x  =  (A  +  BK)x. 


Theorem  1  ('Robust. ness  Form V  Let,  R  >  0  he  fixed.  Suppose  that  (A,  D)  is  symmetrically  stabi¬ 
lizable  wit  hin  A  C  //2",  and  let 

a  :=  sup  max  (A,). 

A6A  !<•'<" 

If  h  G  Sij(|or|)  (resp.  B(\a\)  ),  then  the  nonlinear  control  system  x  =  Ax  +  Bu  +  h(x)  is  stabilizable 
within  B(R)  (resp.  globally  stabilizable)  using  linear  state  feedback. 

Under  the  foregoing  assumption,  this  asserts  the  existence  of  a  feedback  gain  matrix  K  G  lRmxn 
for  which  the  associated  RAS  contains  B(R),  for  each  h  G  R/j(|a|).  No  h  G  Rii(|a|)  results  in  an  RAS 
not  entirely  containing  B(R). 


3.  PLANAR  SYSTEMS 


In  this  section,  we  apply  the  results  of  Section  2  to  planar  systems,  i.e.,  systems  as  in  (1)  with 
A  G  IR2*2  and  b  G  IR2  \  {0}.1  Let 

A  :=  {  (Ai,A2)  G  IR-  x  IR_  :  Aj  ^  A2  and  Ai,A2  <r(A)  }, 

,(AM:=^^b  and  KA,b) 


\b\ 


K 


Also,  define  the  set  A (A,b)  by 


A  (A,  b )  =  {(Aj ,  A2)  G  A  :  Ar  A2  -  u(A,  6)(Aj  +  A2)  +  /i(A,  6)  =  0}, 


(7) 


where  Adj(A)  denotes  the  adjugate  of  A.  The  defining  equation  in  (7)  is  that  of  a  hyperbola.  This 
hyperbola  is  equivalently  characterized  by 


_  u{A,b)\x  -  fi(A,b) 
2  Ar-i ,{A,b) 


(8) 


For  the  next  result,  we  assume  that  the  hyperbola  (8)  is  nondegenerate,  i.e.  it  is  not  a  horizontal 
line.  It  will  be  seen  shortly  that  this  assumption  amounts  to  (A,  b )  being  controllable. 

Theorem  2.  Assume  that  (A,  b)  is  controllable.  Then  the  pair  ( A,b )  is  normally  stabilizable  within 
A(A,6)  if  and  only  if  v(A,b)  <  0.  Furthermore,  given  any  set  of  desired  closed-loop  eigenvalues 
(Ai,A2)  G  A (A,b),  the  corresponding  normally  stabilizing  feedback  gain  is  given  by 


k=[  1  1]  [  (A  J-Ay'b  (A  2I-A)~1b}-1. 


(9) 


Px-oof.  First,  we  show  that  A(A,b)  /  0.  From  Eq.  (8),  we  obtain 

<9Ai  (Ai  -  u(A,  6))2 

1  If  we  let  B  be  a  nonzero  2x2  matrix,  then  it  is  either  nonsingular,  in  which  case  the  stabilization 
problem  becomes  trivial,  or  of  rank  one.  The  latter  case  is  equivalent  to  considering  B  to  be  a  vector  b 
in  IR2. 


Define  the  controllability  matrices  C,  (7a  by  (7  :=  [6  Ab\,  (7a  :=  [ b  Adj(A)6]  and  note  that  det(Ca)  = 
—  det(C')  7^  0.  Then  it  easily  follows  that 

det(ClCa)  =  (det(Ca))2  =  |6|4  (,/(A,fe)  -  v2(A,b))  >  0.  (11) 

and  hence  that  A2  is  a  monotonically  strictly  increasing  function  of  Ai.  A  quick  sketch  of  the  plot  of 
A2  as  a  function  of  Ai  convinces  us  that  A(A,  b)  fl  IR2_  ^  0  precisely  when  u(A,  b )  <  0.  The  sketch  just 
referred  to  is  also  useful  in  finding  pairs  (Ai,  A2)  €  A(A,6).  To  find  such  a  pair,  we  may  simply  pick  a 
value  Ai  in  (i/( A ,  b),  0)  \  cr(A),  and  then  use  Eq.  (8)  to  compute  the  corresponding  value  of  A2. 

Next,  we  show  that  for  (A  1 ,  A2)  €  A(A,  b),  the  vectors 

Vi  =  (A il  —  A)~lb  ,  >  =  1,2 

are  orthogonal.  It  is  equivalent  to  show  that  the  vectors  uq  and  u>2  are  orthogonal,  where 

Wi  =  XA(^i)  Vi  ,  >=1,2 

and  Xa(s)  denotes  the  characteristic  polynomial  of  A.  (Recall  that  A1(  A2  £  <r(A).)  We  have 

wf  w-2  =  b1  (Adj(Ai/  -  A)f  Adj(A2/  -  A)b 

=  bT  (AiAa/  -  Ax  Adj(A)  -  A2AdjT(A)  +  Adjr(A)Adj(A))  b 
=  \bf  AiA2  -  (Ai  +  A2)6TAdj(A)6  +  |Adj(A)6|2  . 

Since  (Aj,  A2)  €  A(A,  b),  it  follows  that  wf  w2  =  0. 

We  now  show  that  iq  and  v2  are  eigenvectors  of  A  +  bk  corresponding  to  Ai  and  A2,  respectively. 
Let  V  =  [vi  n2].  Then  k  =  [1  1]V~ 1  and 

[  ( A  +  bk)vx  (A  +  bk)v3  ]  =  (A  +  bk)V  =  AV  +  bkV 

—  [Av  1  +  b  Av2  +  6]. 

On  the  other  hand,  Avi  +  b  —  (A(A,I  —  A)-1  +  /)  b  =  A,t>,  for  i  =  1, 2.  Therefore 

{A  +  bk)vi  =  A  iVi,  i  =  1,2, 

i.e. ,  a(A  +  bk)  =  { Ai ,  A2}.  Since  the  eigenvectors  vi  and  v2  are  orthogonal,  we  obtain  that  A  +  bk  is  a 
normal  matrix  (in  fact  symmetric  since  Ai  and  A2  are  real). 

□ 

Theorem  2  states  that  if  (A,b)  is  controllable  and  v(A,b)  <  0,  then  every  pair  (Ai,A2)  £  A(A,b) 
may  be  assigned  via  linear  linear  feedback  while  achieving  the  normality  requirement.  It  is  not  necessary 
however  that  a  pair  (A,  6)  be  controllable  for  it  to  be  normally  stabilizable. 

Remarks. 

1.  It  is  easily  shown  that  by  allowing  complex  eigenvalues  in  A  and  A(A,  b),  one  obtains  one  additional 
pair  of  assignable  eigenvalues;  namely  iz(A,6)±  >'|de<(C)|.  Thus  the  set  of  all  distinct  assignable 
eigenvalues  not  in  <r(A)  is  essentially  real. 


2.  The  condition  u(A,b )  <  0  implies  that  det[fc  Adj(A)fc]  yf  0,  hence  (A,  b)  controllable  since 
det[f>  Ab ]  =  — det[A  Adj(A)fr],  in  all  but  the  case  when  b  and  Adj(A)6  are  of  opposite  directions. 

3.  Note  that  it  is  not  possible  to  force  both  eigenvalues  to  be  arbitrarily  large.  Clearly,  this  is  a 
consequence  of  the  normality  requirement.  It  can  be  shown  that  in  order  for  a  pair  (A,B)  to 
be  “arbitrarily”  normally  stabilizable  (normally  stabilizable  with  arbitrarily  negative  assignable 
closed-loop  eigenvalues)  it  is  necessary  and  sufficient  that  rank  ( B )  —  n,  n  being  the  size  of  A. 
The  next  theorem  is  a  direct  consequence  of  Theorems  1  and  2.  Let  R  >  0  be  a  fixed  number. 
Theorem  3:  Assume  that  is(A,b)  <0.  If  h  E  Bh(|^(A,6)|)  (resp.  B(\u(A,b)\)  )  then  there  exists 

a  linear  feedback  u(x)  =  kx  such  that  (the  origin  of)  the  closed-loop  system  x  =  (A  +  bk)x  +  h(x)  is 
asymptotically  stable  within  B(R)  (resp.  globally  asymptotically  stable).  Furthermore,  for  any  desired 
closed-loop  eigenvalues  Ai  e  (  v(A,  b),  -||h||R  )  \  <z(A)  (resp.  (v(A,  b),  -||h|| )  \  cr(A)  )  and  A2  given  by 
(8),  the  feedback  gain  k  is  given  by  (9). 

The  analogue  of  Theorem  1  (Robustness  Form)  in  the  two-dimensional  case  follows  by  taking  a  to 
be  any  number  strictly  between  v(A,  b )  and  0.  The  desired  feedback  k  is  obtained  in  a  manner  identical 
to  that  outlined  in  Theorem  3,  with  -||h||fl  (resp.  -  ||/i||)  replaced  by  a. 

4.  EXAMPLES 

Example  1:  Let  R  ~  1  and  consider  the  system 

•  3  2 

Xl  =  - -X 1  +  x2  +  x{ 

x2  =  x2  +  u  —  x\. 

The  origin  of  the  unforced  system  is  unstable  since  cr(A)  =  {  — §,  1}.  It  is  easily  checked  that  (A,  6)  is 
controllable,  i/(A,  b)  =  -§,  ||/i||fi=1  =  1  <  |  v(A,  6)|  and  ^(A,  b )  =  By  picking  Ai  in  (-§,  )  \ 

<t(A)  =  (— |,  -1),  say  Ai  =  — |,  we  get  from  (8)  that  A2  =  -•y  and  from  (9)  that  k  =  [1  -  2£],  The 
closed-loop  system  is 

3  2 

*1  =  —  2*1  +x2  +  xl 

21  , 
x2-x1-  —X2  -  x2. 

Simulations  of  the  closed-loop  system  for  various  initial  conditions,  shown  in  Fig.  1,  corroborate 
the  fact  that  B(l)  is  contained  in  the  actual  RAS.  Note  that  some  initial  conditions  in  the  immediate 
vicinity  of  B(l)  lead  to  instability  (e.g.,  xo  =  (1.4,0)  and  (1.2, 1.2)). 

Example  2:  We  consider  globally  stabilizing  the  system 

3 

*'i  =  -Tj*!  +  X2 

x2  —  —x\  +  x2  +  u  +  sin(xi). 

Here,  A  and  6  are  as  in  Example  1,  h(x)  =  [0,sin(xi)-xi]T  and  ||/i||  =  1.217.  Proceeding  as  in  Example 
1,  we  obtain  by  choosing  Aj  =  — |  that  the  closed-loop  system  is 

3 

*i  =  —  2*i  +  X2 

21  ... 
x->  =  — “X2  +  sin(xi). 


With  the  Lyapunov  function  V(r)  =  x\  I  x\,  we  line)  Uial.  V {x)  <  -x'f  -  <  0  for  all  x  /  0. 

Therefore,  the  null  .solution  of  the  closed-loop  system  is  globally  asymptotically  stable,  as  predicted  by 
Theorem  3. 


5.  CONCLUSION 

Sufficient  conditions  for  stabilizability  of  nonlinear  systems  with  a  region  of  asymptotic  stability 
containing  a  prescribed  ball  in  IRn  have  been  presented.  Under  a  symmetric  stabilizability  condition  on 
the  system  linearization,  it  was  shown  that  there  is  a  linear  stabilizing  controller,  and  that  the  closed- 
loop  system  stability  is  robust  to  certain  model  perturbations.  Necessary  and  sufficient  conditions  for 
normal  stabilizability  of  a  two-dimensional  linear  time-invariant  system  were  obtained.  These  facilitated 
identification  of  a  closed-form  formula  for  a  stabilizing  feedback  gain  k  which  guarantees  stabilization 
within  a  given  ball.  Necessary  and  sufficient  conditions  for  normal  stabilizability  of  a  general  pair  ( A,B ) 
is  an  issue  that  is  currently  being  considered. 
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